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Abstract 

o ■ 

^vj We investigate a parametric extension of the classical s-dimensional 

Halton sequence, where the bases are special Pisot numbers. In a one- 
Q , ■ dimensional setting the properties of such sequences have already been in- 

^* vestigated by several authors 5, 8, 23, 28 . We use methods from ergodic 

theory to in order to investigate the distribution behavior of multidimen- 

^^ ■ sional versions of such sequences. As a consequence it is shown that the 

Kakutani-Fibonacci transformation is uniquely ergodic. 



Ctf 



1 Introduction 



r*] ' In this article we consider the distribution properties of deterministic point se- 

quences in [0, l) s . We use the following notation: for two points a, b £ [0, l) s we 
write a < b and a < b if the corresponding inequalities hold in each coordinate; 
furthermore, we write [a, b) for the set {x £ [0, l) s : a < x < b}, and we call 
such a set an s-dimensional interval. Moreover we denote by 1/ the indicator 
function of the set I C [0, l) s and by X s the s-dimensional Lebesgue measure, 
for short we write A instead of Ai. Note that vectors will be written in bold 
^J. . fonts and we write for the s-dimensional vector (0, . . . , 0). 

\o ■ 

A sequence (x n )neN of points in [0, l) s is called uniformly distributed modulo 1 

Tt : ( u - d -) if 

m ■ J im at = A s (a,b)) 

for all s-dimensional intervals [a, b) C [0, l) s . A further characterization of 
uniform distribution is due to Weyl 1301 : a sequence (a;„) n6 N of points in [0, l) s 

X 
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is u.d. if and only if for every continuous function / on [0, l) s the relation 

lim £n=l Z( X ") = f /(x)(fx 

holds. Weyl's criterion suggests a numerical integration technique which is usu- 
ally called Quasi Monte Carlo (QMC) integration. However, the theorem gives 
no information on the quality of the estimator. 

The Koksma-Hlawka inequality [18 a states that the integration error of QMC 
integration can be bounded by the product of the variation of / (in the sense 
of Hardy and Krause), denoted by V(f), and the so-called star-discrepancy D^ 
of the point sequence (x n )„ e N, i.e. 



i N r 

-]T/(x„)-/ /(x)d» 



< V(f)D* N (x n ), 



where D* N is defined by 



D* N = D^(xi,...,xjv) = sup 

aG[0,l) s 



ELl(0,a)(xJ_ Aj([())a)) 



N 



In order to minimize the integration error we have to use point sequences with 
small discrepancy. There are several constructions for sequences which achieve a 
star-discrepancy of order (D(N~ 1 (\ogN) s ), so-called low-discrepancy sequences. 
Note that this convergence rate, which is best possible among known sequences, 
is for all s > 1 better than that of the probabilistic error of the standard Monte 
Carlo method, where a sequence of random instead of deterministic points is 
used. QMC integration is successfully applied in several different areas of ap- 
plied mathematics, for example in actuarial or financial mathematics, where 
high-dimensional numerical integrals appear frequently, see e.g. [41 124] . 

In this article we will construct point sequences by a combination of methods 
from uniform distribution theory and dynamical systems. 

Definition 1. Let (X, A, fi) be a probability space. A measurable transformation 
T : X — > X is called ergodic if for every A € A such that T^ 1 (A) = A, either 
fi(A) =0 or n{A) = 1. 

The system [X, A, fi, T) is called a measure theoretical dynamical system, or 
dynamical system, for short. If T is ergodic, the system is called ergodic. 

The link between dynamical systems and uniform distribution is given by the 
following classical result of Birkhoff. 

Theorem 1 (Pointwise Ergodic Theorem). Let (X,A,fi,T) be a dynamical 
system. Then, for every f € £ 1 (X) 

N-l 

lim — V fiT^x) 

3=0 



exists for pL-almost every x G X (here T°x = x). Furthermore if T : X —¥ X is 
ergodic, then for every f G £ (X) we have 

J im m E ^ T3x ) = / fWM*) > 

JV-s-oo iv *■ — ' / v 

j=a JX 

for ^.-almost every x G X . 

Thus if (x„)„ e N is constructed as the orbit of a point x G [0, l) s with respect to 
an ergodic transformation T on [0, l) s , i.e. (x„)„ e N = (T"x) n(E N, we can easily 
see the connection to Weyl's criterion. Nevertheless, the statement of Theorem 
[TJ is only valid /i-almost sure, so we need stronger conditions in order to assure 
that the orbit of a certain (or every) point x G [0, l) s under the transformation 
T is u.d. 

Definition 2. A continuous transformation T : X — ► X on a compact metriz- 
able space X is called uniquely ergodic if there is only one T -invariant Borel 
probability measure on X. The system (X,A,(J.,T) is called uniquely ergodic. 

If the transformation T is uniquely ergodic, then Theorem Q] holds for every 
x G X and so the associated sequence (T n x) ne n is u.d., see e.g. [3S1 Theorem 
6.19]. 

A way to further analyze the properties of a dynamical system is to find an 
isomorphism with a well-known ergodic or uniquely ergodic system. 

Definition 3. Two dynamical systems (JQ, Ai, /J>i,Ti) (i = 1,2) are called iso- 
morphic if there exist two sets A\ G A\ and A^ G A2, with /Ltj(Aj) = 1 [i = 1,2) 
and a bijection (p : A\ — > A2 such that 

T 2 o cp(x) — (p o Ti(x) , Va; G A\ . 

A classical example of a uniquely ergodic system is (Z(,,t) (see for instance, 
[33]), where Zj, is the compact group of 6-adic integers and r : Zf, — > Z^ the 
addition- by-one map (called odometer). We will shortly recall the connection 
between (Z(,,r) and low-discrepancy sequences on [0,1). For an integer b > 2, 
every zgZj has a unique expansion of the form 

j>0 

with digits Zj G {0, 1, . . . , b — 1}. For z G Zf, we define the 6-adic Monna map 
ipb : Zf, — > [0, 1), see also fT3], by 



\i>o / i>o 

The restriction of yt to No is called radical-inverse function in base b and the 
sequence 

(tPb(n))n£N 



is the so-called van der Corput sequence in base b which is a low-discrepancy 
sequence in [0, 1). 

The Monna map is continuous and surjective but not injective. In order to 
make it an isomorphism we only consider the so-called regular representations, 
i.e. representations with infinitely many digits Zj different from 6—1. So the 
Monna map restricted to these regular representations admits an inverse (called 



pseudo-inverse) 


tf: 


[0, 


1)- 


— > Zj,, 


defined by 










tf 


\j>0 


*b-*A = 


3>0 


where J2j>o z j^ 


-j'-i 


is 


a 6-adic rational in [0. 


,!)• 



Moreover tpi, is measure preserving from Zj, onto [0, 1) and it transports the 
normalized Haar measure on Zj, to the Lebesgue measure on [0, 1). Hence by 
the unique ergodicity of r it follows that the sequence (t"z)„>o is uniformly 
distributed in Zj, for all z £ Z&, in particular for z = 0. Thus the van der 
Corput sequence (^b(r™0)) ra6 N in base b is uniformly distributed modulo 1. 

In order to construct multidimensional sequences we need a criterion to ensure 
that the Cartesian product of several ergodic systems is again ergodic, see e.g. 

D3I- 

Theorem 2. Let % — (Xj, Ai, jjl%, Tj), i = 1, . . . , s, be uniquely ergodic dynam- 
ical systems. Then the dynamical system 7i X ... X T s is uniquely ergodic if and 
only if for all i,j £ {1, . . . , s}, i ^ j, the discrete parts of the spectra of 7* and 
Tj intersect only at 1. 

In 13 , the spectrum of (Z{,,r) is given explicitly. Furthermore the authors 
show that if b\ , bi denote positive integers and t^ denotes the addition- by- 
one on Zb i , then the dynamical systems Zjj and Z& 2 are spectrally disjoint 
if and only if b\ and 62 are coprime. This is exactly the condition proved 
by Halton in |17) in order to obtain a low-discrepancy sequence in [0, l) s by 
combining coordinatewise van der Corput sequences. The resulting sequence 
(0b(n)) ? ieN = (<t>bi(ri), ■ ■ ■ j 4 > b d (n)) n & : i is called b-adic Halton sequence, where 
b is an s-dimensional vector of pairwise coprime integers bi, i = 1, . . . , d. 

The aim of the present article is to extend the above idea to point sequences 
with irrational bases. Such sequences in the unit interval were investigated by 
several authors. Barat and Grabner |5_ consider the so-called /3-adic van der 
Corput sequence (<^(n)) rag N on [0,1) and similar constructions. They prove 
that (</>/3(n))neN is low-discrepancy, where j3 is the characteristic root of special 
linear recurrences. Ninomiya [23] considers the discrepancy of point sequences 
on [0, 1) for a slightly greater class of irrational bases j3. In [20] the underlying 
construction was extended to piecewise linear maps. Furthermore Steiner [27] 
considers so called bounded remainder sets in this setting. In a second article 
Steiner [28] considers van der Corput sequences on abstract numeration systems 
and gives conditions under which they are low discrepancy. Discrepancy bounds 



for a higher dimensional extension of [53; are given in [33]. Note that the con- 
struction in |22| is different from that in the present article. 

Carbone [5] and Drmota and Infusino [TU] investigate the discrepancy of point 
sequences generated by the so-called Kakutani splitting procedure. Carbone 
completely characterises the growth order of the discrepancy for a two para- 
metric subfamily, the so-called LS-sequences. Moreover, Aistleitner, Hofer and 
Ziegler [T] give conditions under which an s-dimensional vector of LS-sequences 
is not u.d. in [0, l) s . 

The remainder of the article is structured as follows: in the next section we 
formulate a characterization of uniquely ergodic systems which are constructed 
as Cartesian products of odometers on numeration systems related to linear 
recurrences. In the third section we give conditions under which a construction 
like that of Halton produces an u.d. sequence on [0, l) s with respect to irrational 
bases f3\ , . . . , f} B . Furthermore we present a parametric class of sequences which 
satisfies these conditions. Finally we prove that the ergodic Kakutani-Fibonacci- 
transformation, presented in [5], is in fact uniquely ergodic. 

2 General G-odometers 

In this section we consider odometers on numeration systems, which are related 
to linear recurrences. For a detailed discussion of such number systems we refer 
to [TTJ [TH [T21 [H] . We first recall some basic results and definitions and then 
investigate Cartesian products of odometers. 

Definition 4. Let (G„)„>o be an increasing sequence of positive integers with 
Go = 1 • Then every positive integer can be expanded in the following way 

oo 

VneN, n = Y,£kG k7 (1) 

fc=0 

where e k £ {0, . . . , [Gk+i/Gk\} and [x\ denotes the integral part of x, that is the 
greatest integer less than or equal to x £ R. This expansion (called G-expansion) 
is uniquely determined and finite, provided that for every K 

A'-l 

J2 £fcG fc < G K . (2) 

fc=0 

The digits £k can be computed by the greedy algorithm (see for instance [TTj ) 
and the sequence G = (G„)„>o is called numeration system. 

We denote by Kq the subset of sequences that verify the property ^ and the 
elements in K.q are called G-admissible. In order to extend the addition-by-one 
map defined on N to ICq we introduce KP G C Kg 

1C g = Ixe)C g :3M x ,Vj>M x ^e fc G fc < G j+1 - 1 1 . (3) 

I k=o J 



If we denote by x(j) — J2k=o £ kGk, then we set 

t{x) = (e (x(i) + 1) . ..£j(x(J) + l))e j+1 e ]+2 ... , (4) 

for every x 6 KP G and j > M x . This definition does not depend on the 
choice of j > M x and can be easily extended to sequences x in Kg \ K G by 
t(x) — = (0°°). In this way the transformation r is defined on ICq and it is 
called G-odometer. We refer to [14] for a complete survey on odometers related 
to general numeration systems. 

In this article we consider only numeration systems where the base sequence is 

a linear recurrence. Let Go = 1 and Gu = aoGk-i H h <2fc_iGo + 1 for k < d. 

Then G n for n > d is determined by a recurrence of order d > 1, i.e. 

G n +d — aoG n+d -i H hfld-iG„ 77, > . (5) 

The solution of the characteristic equation of the numeration system G 

x d = aox*- 1 + . . . + a d - x . (6) 

plays a central role. We will be mainly interested in numeration systems where 
the solution of ^ is a Pisot number j3. Note that f3 is always a Pisot number if 

a > ... > a d -i > 1, (7) 

see Theorem 2]. By [2S] we get that in this case the so-called Parry's (3- 
expansion of f3 is finite, i.e. 

a i ai i i flrf - 1 ?q\ 

P = a + J + ... + WT , (8) 

where ao = |_/5 J . At the end of the last section we will also consider numeration 
systems where ([7|) does not hold. 

For numeration systems where the characteristic root /3 is a Pisot number which 
satisfies (|5|), we have that a finite sum Xlfe=o £fe ^ fe ' s ^ ne expansion of some 
integer if and only if the digits Sk of the G-expansion satisfy 

(e^efe-i,...^^ 00 ) < (a ,ai,...,a d _i)°° , (9) 

for every k and < denoting the lexicographic order (see [25]). Representations 
(ek, ■ ■ ■ , So) verifying this condition are called admissible representations and so 
they belong to /Co- 
in [2] Theorem 5] , the authors show that the odometer on an admissible numer- 
ation system G is uniquely ergodic and that the corresponding unique invariant 
measure fj, is given by 

lx{Z) = (10) 

Fk^- 1 + (F K+1 - a F K )p d - 2 + ... + (FK+d-j - a F K+d - 2 - ... - a d - 2 F K ) 

P K (P d -i+P d - 2 + ... + l) 

where Fk :— #{n < Gk '■ n G Z} and Z is the cylinder with fixed digits 
eo, . . . , £k—i- Note that the formula in [141 Theorem 5] included a misprint and 



was stated in corrected form in [5]. 

In the sequel we want to apply Theorem (5J thus we need information on the 
spectrum of the G-odometer. We introduce the following two hypotheses: 

Hypothesis 1 (Grabner, Tichy and Liardet [14] U There exists an integer b > 
such that for all k and 



^ 



i=0 j=k+b+2 

the addition of G m to N , where m > k + b + 2, does not change the digits 
eo, . . . , Cfe, in the greedy representation i.e. 

k oo 

Hypothesis 2 (Frougny and Solomyak [T2]). TTie solution j3 of equation ^ is 
a Pisot number such that all numbers of the set Z[/3 _1 ] have finite (3- expansions. 

In |14j the authors remark that the Multinacci sequence, i.e. oo = • . • = dd—i = 
1, fulfills Hypothesis [TJ Several authors worked on algebraic characterizations 
of Pisot numbers j3 which satisfy Hypothesis [2] Frougny and Solomyak show 
that ^ implies Hypothesis [5] and they give a full characterization of all Pisot 
numbers of degree two with this property. Furthermore Hollander |19] states 
another sufficient condition for Hypothesis [5] and Akiyama [5] characterizes all 
Pisot units of degree three satisfying Hypothesis [5] Further progress was also 
made by Akiyama et al. [3] who prove Hypothesis [2] for a large class of Pisot 
numbers of degree three by using the theory of shift radix systems. Nevertheless 
there exists no complete algebraic characterization for Pisot numbers satisfying 
Hypothesis [5] of degree greater than two. Note that both hypotheses can be 
satisfied by the same numeration system but, to the best of the authors knowl- 
edge, it is unknown if the two hypotheses are equivalent, see [14]. 

Grabner, Liardet and Tichy [T2J Theorem 6] and Solomyak [^Hl Theorem 4.1] 
show that the odometer on the base system G has purely discrete spectrum 
provided that one of the above Hypotheses holds. Furthermore we obtain in 
both cases that the set of eigenvalues of the transformation is given by 

r := {z e C : Km z G " = 1}. (11) 

n— >oo 

Theorem 3. Let G , . . . , G s be numeration systems given by ([5]). Assume that 
the coefficients of the linear recurrences are given as ah — bi, i = 1, . . . , s, j — 
0, . . . , (di — 1), with pairwise coprime, positive integers bi,i = l,...,s. Fur- 
thermore let %■ ^ Q, for all l,k e N, where (3i,...,f3 s are the roots of the 

characteristic equations ©. Then the dynamical system which is constructed 
as the s-dimensional Cartesian product of the corresponding odometers, i.e. 
((/Cgi,ti) x ... x (ICg^Ts)), is uniquely ergodic. 



Proof. It follows by [26, Main Theorem] that the G % fulfill Hypothesis [5] and thus 
the components of the s-dimensional dynamical system are uniquely ergodic. 
Furthermore we obtain that their spectrum is given by (|11[) . By Theorem [51 we 
derive that the Cartesian product is uniquely ergodic if and only if T t n Tj = 1 
for all 1 < i < j < d. As noted in [14] , we have the following connection between 
Pi and the corresponding sequence G l n , 

Hm^=C„ (12) 

where the constant C{ can be computed by residue calculus. Now consider a 
fixed I £ N and 

exp ( 27ri^f J a exp faiC^' 1 ) 

V Pi J 

« exp (27r?G^_;) , 
and thus 



. G^ 



lim^exp ( 27ri^- ) = lim^exp (2iriG l n _ l ) = 1, 



where G^ is given in (TT21 . Furthermore, it is easy to see that for every k £ N 
there exists a uq with 6^ | G„ for all n > tiq and there exist no b',n' £ N with 
gcd(6', 6j) = 1 such that 6' | G n for all n > n' . By simple considerations we get 
that Tj can be written as 

exp (271-z-^-r | : m,Z,c, e NU{0} 

An application of Theorem [5] completes the proof. D 



3 Uniform distribution of the /3-adic Halton se- 
quence 

The purpose of this section is to formulate conditions on the odometers {JCqi ,t±), 
. . . , (/Cg s 7 T s) such that their product dynamical system is uniquely ergodic and 
the Monna map transports the measure fj,i X . . . X fj, s , where the (ii are given 
by (fTO)) . to the Lebesgue measure on [0, l) s . Under such assumptions we show 
the resulting s-dimensional, /?-adic Halton sequence to be u.d. in [0, l) s . 

First we extend the definition of the Monna map to irrational bases /3 > 1. Let 

n = ^Z € i G J 
i>o 

be the G-expansion of an integer n. We define the /3-adic Monna map (f>p : K.q —> 

R+ as 

\j>0 J j>0 



We call 

j>0 

the /3-expansion of x. Furthermore, as in the first section, we define the radical 
inverse function as restriction of (f>p on /C^ and define the pseudo-inverse <pt 
similarly. In this context we define the /3-adic Halton sequence as 4'f3( n ))neN — 
(4'j3 1 (n), . . . ,<j)/3 s (n)) n£ ^, where f3 = (/3i, . . . ,/3 s ) and the /3i are solution of the 
corresponding characteristic equations. 

Note that even if a Pisot number j3 is chosen as the solution of ©, it is not 
sure that the image of /Co under <f>p is a subset of [0, 1) or dense in [0, 1). The 
following lemma gives a characterization of numeration systems for which this 
is true. 

Lemma 1. Let a = (oo, . . • ,a,d—i), let the integers ao, . . . ,ad—i > be the 

coefficients defining the numeration system G and assume that the corresponding 
characteristic root fj satisfies ([B]). Then <fip(N) C [0,1) and <f>p(N) (jL [0,x) for 
all < X < 1 if and only if a can be written either as 

a=(oo,...,oo), (14) 

a = (a ,a - 1,. . . ,a - l,a ), (15) 

a = (a ,...,a ,a + 1) (16) 

or 

a=(a',...,a',a"), (17) 

where a$ > 0, a', a" are of equal length and are of the form 

a' = (a , . . . , a , a - 1), a" = (a , . . . , a ) or 

a' = (a a ,a Q - 1, . . . , a - 1), a" = (a ,a Q - 1, . . . , a - 1, a ) 

Proof. It follows by @ that 



[3 ••• 0" 



(18) 



Furthermore we know that for all admissible representations of an integer n we 
have 

(sfejEfc— l) ■ • • , £ o,0°°) < (ao,ai, . . . ,ad_i)°° , 

for every k and < denoting the lexicographic order. If we have ao ^ maxo<i<d(ai), 
by ([T5)l we obtain 

CO 

max(A«(n) < > — - < 1, 

i=l H 

where equality holds for both inequations only if a is of the form (TT6")) . 

Assume now that ao = maxo<i<d(ai) and there exist an < j < d such that 
a,j < ao — 1 and k > is the maximal integer with a, = ao for all i < k. Then 
the representation < k < d — 1 and (ao, . . . , afc, (ao — 1, ao, . . . , afc_i)°°) is 
admissible. But by flT8"|) we get that the image of the /3-adic Monna map of this 
representation is strictly greater than 1. By a similar argument can be applied 



in the case a c i-i 7^ ao- 

Note that we have excluded every case which is not of the form (ao, fflj., . . . , ad-2, ao), 
where aj 6 {ao — l,ao} for 1 < « < d — 2. By (ITS)) and the construc- 
tion of the van der Corput sequence we obtain that the image of the /3-adic 
Monna map of (ao, . . . , a&-i-, a<i-\ — 1)°° is 1. Thus it is sufficient to show that 
(ao, . . . , ad-2, Q>d-\ ~ 1)°° is admissible and its image under the Monna map is 
maximal. This is clear if a is of the form (fTl)| - ([T7|) . 

Finally assume that a = (ao, 01, . . . , ad-2, ao), where a* £ {ao — 1, ao} for 1 < 
i < d — 2, a is not included in one of these cases and let k be defined as 
above. Then < k < d — 1 and (ao, . . . , a*,, (ao — 1, ao, . . . , afc_i)°°) is maximal 
and admissible. But this is only equal to (ao, . . . , ad-2, a>d—x ~ 1)°° when (|17[) 
holds. D 

Remark 1. Note that (|17|l is another way to represent number systems defined 
by a" satisfying (1141) or (|15[) . For example the numeration system defined by 
a* = (1, 0, 1, 1) is the same as the Fibonacci numeration system where a = (1, 1). 
Furthermore if a is of the form (|16p we can rewrite it as the classical ao-adic 
numeration system. Hence in the sequel we will only be interested in numeration 
systems which fulfill flU) or ((To)) . 



Lemma 2. Let G be a numeration system of the form ([5]). assume that the 
coefficients of the linear recurrence are given by a,j = a, j = 0, . . . , (d— 1), for a 

positive integer a and let [3 denote the corresponding characteristic root. Then 
n(Z) = \(4>p(Z)) for every cylinder set Z. 

Proof. Let the cylinder set Z be defined by the fixed digits eo, . . . , efc— l- Assume 
first that e fc _i < a, then F k+r = (a + l) r for < r < d. Thus, by [[ID]), we 
obtain that 

M(£) = /r fe - 

Consider the /3-adic Monna map at n £ N, i.e. 

i=0 ^ 

If efe_i < a we can easily see that 4>p{Z) is dense in 

E e» v^ e^ (efc-i +1) \ 

L»=0 ^ j=0 ^ r / 

and that <f)p(x') £ I if x' £ Z. Thus 4>p(Z) is A-measurable and A(<^(Z)) = 
A(7) = /?- fc . 

Assume now that ^ is defined by the fixed digits eo, . . . , ek-2 and e^-i = a. By 
the above argument we derive that a cylinder with fixed digits eo, . . . , £k-2 has 
measure /3~( fc_1 ) and every cylinder with digits eo, . . . , e^-i has measure (3~ h . 
Thus 

H(Z) - i3- {k ~ l) -(a- f)r k - 



10 



Now we consider <j)p(Z), hence 

"fe-l fe-3 



MZ) 



and thus \((ftp(Z)) = fJ,(Z) 



V" e * , ( e ^-2 + !) 



,i=0 M i=0 K K 



Let 2 < h < min(/c, d — 1) and consider a cylinder set Z with fixed digits 
eo, • • • , tk-h-i < a an d £k-i = a iov I = 1, . . . , h. Then, as above we get that the 
cylinder with fixed digits eo, . . . , €k~h-i has measure (3~( k ~ h > and every cylinder 
with digits eo, • . • , ek-h+i has measure fj-( k - h + 2 ) _ Thus we get that 

n(Z) = /r < fc -" +1 ) - (a - i)/r ( fe -"+ 2 > 

Considering (ftp(Z) we have 

'k—h+l k—h 



h(z) 



2-~i Ri+1 ' 2-~i 



£i . (Cfc-h-1 + 1) 



Ri+1 ' L^ Ri+1 ak-h+2 

i=0 H i=0 H H 



and thus \{<j>p{Z)) = \x{Z). D 

Remark 2. As mentioned in the previous section, a result of Frougny and 
Solomyak ] 12j Lemma 3] implies that the dominant root of 

2 \ i 

x — a$x — ox, ao,ax>l, 

is a Pisot number if and only if ao > a\. By Lemma U we know that the 
image of K.q under (ftp is not a subset of [0, 1), when ao > ax- Thus Lemma [2] 
characterizes all van der Corput-type constructions when d = 2. 

Theorem 4. Let G , . . . ,G S be numeration systems as in Theorem^ and let 
j3i, . . . , p s denote the roots of the corresponding characteristic equations. Then 
the s- dimensional, (3-adic Halton sequence (</>/3(ro)) n6 N is u.d. in [0, l) s . 

Proof. By Lemma [2] and the definition of the Monna map we obtain an isometry 
between the dynamical systems {{K.qi ,t\) x ... x {1Cg s ,t 8 )) and (([0, 1), T\) x 
...x ([0,1), T s )) where 

Ti : [0, 1) -> [0, 1), Ti{x) := (ftp z o n o 0+ (x). 

Let Tx = (TxXx, . . . ,T s x s ) for x = (xi, ...,x s ) £ [0, l) s . Hence by Theorem 
[U (T™x)„ e N is u.d. in [0, l) s for all x e [0,l) s . In particular (<f>0(n)) n &n = 

(T"0) neN is u.d. D 

Remark 3. Note that the classical 6-adic Halton sequence with pairwise coprime, 
integer bases bi, . . . ,b s > 2, is included in Theorem @] 

Theorem 5. Let the numeration system G be defined by the coefficients (ao, a\, 02) 
(1, 0, 1) and let j3 be its characteristic root. Then fi(Z) = \{(ftp(Z)) for all cylin- 
der sets Z. Thus T(x) — (ftp o r o <ftt(x) is uniquely ergodic and (T n x)„ S N is 
u.d. for all x in [0, 1). Furthermore the spectrum of T is given by 

r= |cxp(27rz-|) : m,/,ceNU{0}|. (19) 
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Proof. Note that j3 is a Pisot number and equation ((5} holds since [_/3J = 1 = oq. 
Hypothesis [1] was proved for this case in [7J Theorem 4] . The proof that Hy- 
pothesis [5] is fulfilled can be found in [21 Theorem 3] . Equation (|19[) follows by 
the proof of Theorem [3] 

Now we have to prove that (j>p transports the measure [i to the Lebesgue measure 
on [0,1). First we assume k > 3. Let the cylinder Z be defined by the fixed digits 
eo,..-,efc_i. We consider four different cases, first e k s = e k - 2 = Cfc-i = 0. 
Then F k = l,F k+1 = 2, F k+2 = 3 and we get by (HHD that 

M (Z) = /T fc - 

Furthermore by the same argument as in the first part of the proof of Theorem 
3] we obtain 



Mz) 



fc-l fe-2 , 1 ^ 



fli+l '.2^ fli+l 
,i=0 H i=0 ^ 



/3 fe 



and thus \{<fip(Z)) = (3~ k . 

Now let e fc _ 3 = l,e fc _ 2 = £fe-i = 0. Hence F k = l,F k+ i = 2,F k+2 = 3 and 

fi(Z) = p- k . We have 



MZ) 



'fc-1 fc-1 oo N 

2_^ fli+l ' 2_^ fli+l + * 2-^p 

i=0 i=0 i=0 / 

'fe-l fe-l \ 

Z^ fli+1 'Z^ fli+1 ^^ I 

i=0 H i=0 P / 



thus we have again \((j>p(Z)) = /3 fc . Now assume efc_2 = ljCfe-i = 0. Hence 
F fc = l,F fc+1 =l,F fe+2 = 2and 



/?" 



1 



^ = ^Vi + /3 -x + l- 



Similarly as above we get 
rfe-i 



MZ) 



fe-i 



E^E^+/^ (fc+1) E/^ +1) 

i=0 P i=0 ' i=0 



'fe-l fe-l 



V - V - ■+- /?-( fe + 1 ) 

j=0 H i=0 M 



thus A(^(Z)) = /T( fc+1 ). Now we obtain 



P' 



-(fe+i) _ ^-fe 



/3" 



1 



p-t + p-i + l 

<* /T 3 + P~ 2 + P~ l = P~ 2 + 1 



which holds by ([5]). 

In the last case we assume e k -i = 1, thus F k = F k+ \ = F k+2 = 1 and 



Kz) = /r fe 



l 

,5- 2 +/3- 1 + 1' 
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As above we get \{4>p{Z)) — (3 ( fc + 2 ) and the result follows by 



p-(k+2) = p-k L 



p-t + p-i + l 

^ /r 4 + /r 3 + /r 2 = 1 

«*■ P~ 3 + /T 1 = 1. 
The cases, where k < 3, follow by the same arguments. □ 

Remark 4. As a consequence of Theorem [2] we can construct two-dimensional 
u.d. sequences (</> / 3 1 (n),^ ( g 2 (n)) T i 6 N, where fti is the characteristic root in Theo- 
rem[Sl (4>p 2 (n)) n( =N is the characteristic root of a numeration system in Theorem 

0]and ~t 4. Q for all integers k,l > 0. In this way we can construct a new class 

of multidimensional u.d. sequences. 

Note that Theorem[5]extends the examples given in Proposition 13,14], where 

the authors consider G-additive functions which lead to u.d. point sequences in 

the unit interval. 

Furthermore, it is possible to show that the one-dimensional point sequence in 

the previous theorem is a low-discrepancy sequence by mimicking the proof for 

the 6-adic van der Corput sequence, see e.g. [2T1I51I8]. 

In [5] , the authors present the so-called Kakutani-Fibonacci-transformation for 
which they show that it is an ergodic transformation on the unit interval and 
the orbit of is exactly the LS-sequence with parameters L — S = 1 . With our 
approach we can show that this transformation is in fact uniquely ergodic, i.e. 
the orbit of x under the transformation is u.d. for every x G [0, 1). 

Theorem 6. The Kakutani- Fibonacci-transformation is uniquely ergodic. 

Proof. In PP, the authors use a van der Corput-type construction to compute 
the points of the LS-sequence with general parameters L, S. One can easily see 
that this construction is equivalent to the construction presented in the present 
paper, when L = S = 1, G„+2 = G n +i + G n , Go = l,Gi =2 and j3 is the 
golden ratio. Thus, by Theorem |4j we get a uniquely ergodic transformation 
T: [0, 1) — > [0, 1) given by T(x) := <j) p o t o cjjt (x) . Similar as for the classical van 
Neumann-Kakutani map, see |13j . we observe that T is exactly the piecewise 
translation map given in [5]. □ 

Remark 5. Note that, since T is uniquely ergodic, the orbit of x is u.d. for every 
x € [0, 1). This fact is used in so-called randomized QMC techniques, where the 
starting point A is a uniformly distributed random variable in [0, 1). This idea 
can of course be extended to the s-dimensional case. For more information on 
this topic and applications in financial mathematics see e.g. |24j . 

References 

[1] C. Aistleitner, M. Hofer, and V. Ziegler. On the uniform distribution mod- 
ulo 1 of multidimensional ls-sequences. to appear: Annali di Mathematica 
Pura et Applicata, 2013. 



13 



[2] S. Akiyama. Cubic pisot numbers with finite beta expansions. In F. Halter- 
Koche and R. Tichy (eds.): Algebraic Number Theory and Diophantine 
Analysis, page 11. 

[3] S. Akiyama, H. Brunotte, A. Petho, and J. M. Thuswaldner. Generalized 
radix representations and dynamical systems ii. Acta Arithmetica, 121. 

[4] S. Asmussen and H. Albrccher. Ruin Probabilities, 2nd eds. World Scien- 
tific, Singapore, 2010. 

[5] G. Barat and P. Grabner. Distribution properties of g-additive functions. 
Journal of Number Theory, 60:103-123, 1996. 

[6] A. Brauer. On algebraic equations with all but one root in the interior of 
the unit circle. Mathematische Nachrichten, 4:250-257, 1951. 

[7] H. Bruin, G. Keller, and M. Pierre. Adding machines and wild attractors. 
Ergodic Theory and Dynamical Systems, 17(6):1267-1287, 1997. 

[8] I. Carbone. Discrepancy of LS-sequences of partitions and points. Ann. 
Mat. Pura Appl. (4), 191(4):819-844, 2012. 

[9] I. Carbone, M. R. Iaco, and A. Volcic. A dynamical system approach to the 
Kakutani-Fibonacci sequence, to appear: Ergodic Theory and Dynamical 
Systems, 2012. 

[10] M. Drmota and M. Infusino. On the discrepancy of some generalized Kaku- 
tani's sequences of partitions. Uniform Distribution Theory, 7:75-104,2012. 

[11] A. S. Fraenkel. Systems of numeration. Amer. Math. Monthly, 92(2):105- 
114, 1985. 

[12] C. Frougny and B. Solomyak. Finite beta-expansions. Ergodic Theory and 
Dynamical Systems, 12:713-723, 1992. 

[13] P. Grabner, P. Hcllckalek, and P. Liardet. The dynamical point of view of 
low-discrepancy sequences. Uniform Distribution Theory, 7(1):11— 70, 2012. 

[14] P. Grabner, P. Liardet, and R. Tichy. Odometers and systems of numera- 
tion. Acta Arithmetica, 70:103-123, 1995. 

[15] P. Grabner and R. Tichy. Contributions to digit expansions with respect to 
linear recurring sequences. Journal of Number Theory, 36:160-169, 1990. 

[16] P. Grabner and R. Tichy. a-expansions, linear recurrences, and the sum- 
of-digits function. Manuscripta Mathematica, 70:311-324, 1991. 

[17] J. H. Halton. On the efficiency of certain quasi-random sequences of points 
in evaluating multi-dimensional integrals. Numer. Math., 2:84-90, 1960. 

[18] E. Hlawka. Funktionen von beschrankter Variation in der Theorie der 
Glcichverteilung. Annali di Mathematica Pura et Applicata, (54):325-333, 
1961. 



14 



[19] M. Hollander. Linear numeration systems, finite beta- expansions, and dis- 
crete spectrum of substitution dynamical systems. PhD Thesis, Washington 
University, 1996. 

[20] Y. Ichikawa and M. Mori. Discrepancy of can der Corput sequences gen- 
erated by piecewise linear transformations. Monte Carlo Methods and Ap- 
plications, 10(2):107-116, 2004. 

[21] L. Kuipers and H. Niederreiter. Uniform distribution of sequences. Wiley- 
Interscience [John Wiley & Sons], New York, 1974. 

[22] M. Mori and M.Mori. Dynamical systems generated by algebraic method 
and low discrepancy sequences. Monte Carlo Methods and Applications, 
18:327-351, 2012. 

[23] S. Ninomiya. Constructing a new class of low-discrepancy sequences by 
using the /3-adic transformation. Math. Comput. Simulation, 47(2). 

[24] G. Okten, B. Tufhn, and V. Burago. A central limit theorem and im- 
proved error bounds for a hybrid-Monte Carlo sequence with applications 
in computational finance. Journal of Complexity, 22(4):435-458, 2006. 

[25] W. Parry. On the /3-expansions of real numbers. Acta Math. Acad. Sci. 
Hungar., 11:401-416, 1960. 

[26] B. Solomyak. Substitutions, adic transformations and beta-expansions, in: 
Contemporary Mathematics (Proc. of the Conference on Symbolic Dynam- 
ics, Yale, 1991), 135:361-372, 1992. 

[27] W. Stcincr. Regularities of the distribution of /3-adic van der corput se- 
quences. Monatshefte Mathematik, 149(1):67-81, 2006. 

[28] W. Steiner. Regularities of the distribution of abstract van der corput 
sequences. Uniform Distribution Theory, 4(2):81-100, 2009. 

[29] P. Walters. An introduction to ergodic theory, volume 79 of Graduate Texts 
in Mathematics. Springer- Vcrlag, New York, 1982. 

[30] H. Weyl. Uber die Gleichverteilung von Zahlen mod. Eins. Mathematische 
Annalen, 77:313-352, 1916. 



15 



